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. -^ ■ Abstract: The usual procedures to extract effective refractive indices for 

*vq \ negative-index metamaterials are complicated by branching ambiguities in 

the inverse cosine function. The existing methods to eliminate ambiguities 
either involve calculations with varying geometries which are inherently 
'^ ■ flawed, as metamaterial properties depend strongly on geometry, or are 

r~| \ more subjective as one must inspect and select branches to yield effec- 

tive parameters as a continuous function of frequency. We propose that 
g_) \ the Kramers-Kronig relations, which hold for negative-index materials, 

naturally guide the selection of the proper branches, and in fact predict 
negative refractive index where other extraction procedures alone might 
C^ ■ not. The experimental realization of high quality negative-index materials 

requires a robust and reliable index retrieval procedure; the combined 
extraction/Kramers-Kronig retrieval can be used to design optimal metama- 
terials in a simple, systemmatic, and general manner. 
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1. Introduction 

Artificial materials known as metamaterials can exhibit a range of unique electromagnetic prop- 
erties that are not observed in naturally occuring materials. Obtaining a negative index of refrac- 
tion is particularly desirable, as negative-index materials (NIM) have many potential applica- 
tions including imaging beyond the diffraction limit and optical cloaking. Although NIM were 
characterized by Veselago HI many years ago, it is only recently that we have seen a sudden 
burst in interest in NIM design and applications. Pendry ||2l showed that perfect subwavelength 
imaging or superlensing would be possible with a NIM of refractive index A^ = — I + 0/. The 
exciting idea here is that imaging would be limited not by the wavelength of light, but by the 
quality of the lense. Double negative materials and NIM were shortly thereafter experimentally 
realized in the work of Smith and coworkers ID ID. In addition, the functional range of NIM 
has been driven to visible frequencies by Shalaev by incorporating the visible-range electro- 
magnetic resonances of metallic nanoparticles into NIM design Q. In practice, sufficient con- 
ditions to yield a negative refractive index come in the form of Re(£) < 0, Re(/i) < 0, where £ 
and jj, represent the relative permittivity and permeability of the material, respectively. Materi- 
als having simultaneously negative permeability and permittivity are known as double negative 
materials. Note, however, that these conditions are sufficient but not necessary to achieve a 



negative index. Necessary conditions require only that Re(e)Im(/i) +Im(e)Re(/i) < ||6l. 

Several classes of NIM have been realized experimentally with a few standard designs emerg- 
ing that yield negative indices in a wide range of frequencies. Photonic crystals QUI ID have 
been shown to exhibit a negative index in the near-IR range, but their utility as imaging devices 
is limited by the internal structure of the crystal, which has dimensions on the order of the 
wavelength of light QOlIIIl. Split-ring resonators (SRR) Ell [H [H [Bl [l6l [17] can exhibit 
negative index at telecommunication wavelengths. In general, the SRR provides a negative per- 
meability that, when combined with the negative permitivity of a metallic component, results 
in a negative real component of the refractive index. The operating frequencies of the SRR can 
be shifted toward the visible rangle by miniaturization. The first class of NIM to reach the vis- 
ible spectrum is perhaps the simplest and in that sense the most versatile. A pair of nanowires 
separated by some dielectric has been shown theoretically and experimentally give a negative 
refractive index in the near IR and visible range ||5l[Ill[l9l|20l. There also exist S-shaped varia- 
tions on the wire-pair structure in the spirit of the SRR known as S-shaped resonators 12111221 . 
The behavior of double S-shaped resonators can be mimicked by a single pair of asymmetric 
parallel nanowires 1231 . which bear strong resemblence to the simple metal-dielectric-metal 
sandwich described above. 

Clearly there exist many possibilities for high quality NIM, and characterizing the refractive 
index theoretically can aid the design of such materials, circumventing the need for expensive 
experimental optimizations of metamaterial geometeries. Refractive index retreival methods 
based on scattering- or S-parameters lT6l l24l l25l l26ll are the most direct and robust methods 
for relating refractive index to theoretically and experimentally tractable parameters such as 
complex transmission and reflection coefficients. The extraction procedures work reasonably 
well in most cases, but, occasionally, severe branching ambiguities muddle the procedure and 
introduce an element of subjectivity into the index assignments. In this paper, we suggest a 
simple and powerful method for automatically selecting the correct branches in the extraction 
procedure by invoking the concept of causality. The well-known Kramers-Kronig relations au- 
tomatically guide the extraction procedure and in some cases allow us to predict a structure to 
be a NIM that would otherwise be difficult to characterize. 

2. Theory 

S-parameter retrieval procedures begin with the assumption that the metamaterial may be ap- 
proximated by a homogeneous slab with a well-defined thickness, D. The complex transmission 
and reflection coefficients are given by Q 
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where 5 is the distance from the center of the slab to the evaluation planes above and below 
the sample. The reflected and transmitted fields are normalized by the incident light, £, at the 
reflection evalulation plane. The refractive index for a homogeneous NIM in vacuum (or on a 
substrate with refractive index equal to unity) is given by 

N = ± — cos"' + . (3) 

kD \ 2t J kD ^ ^ 

The sign in Eq. Q is uniquely determined by requiring that Im(N) be positive. However, the 
real portion of the inverse cosine function is prone to branching ambiguities, where seemingly 



valid solutions exist for all integers, m. The S-parameter retrieval is based in the assumption 
that the metamaterial may be treated as a homogeneous material, and as such, the macroscopic 
refractive index should be independent of the thickness of the NIM, D. To maximize the separa- 
tion of the branches, the extraction procedure should be performed on as thin a slab as possible, 
and even then, one should perform more than one analysis with differing values for D. The 
true branches should reveal themselves in their insensitivty to the choice of D. Unfortunately, 
metamaterial properties are very sensitive to geometric parameters, and it is virtually impos- 
sible to alter D in such a way as to maintain a constant refractive index. Choosing branches 
becomes a more subjective process wherein one selects parameters as a continuous function of 
frequency lfT6ll24l . We will show that this process is still flawed, as the solution can be skewed 
by specious discontinuities in the refractive index as a function of frequency. 

The ambiguities in the extraction procedure are ultimately rooted in (and may ultimately be 
circumvented by) the fact that Eq. ^ is insufficient to guarantee that the obtained refractive 
indices are causal; that is, do the real and imaginary components reflect the analytical properties 
of e and jj. ? The well-known Kramers-Kronig (KK) relations 1281 |29l guarantee that knowledge 
of the imaginary component of the refractive index over all frequencies implies knowledge of 
the real component by the integration 

«(co)-„. = lpr^;^dco'. (4) 
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Similarly, knowledge of the real component over all frequencies can be mapped to the imagi- 
nary component: 
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Regardless of the extent to which branching issues complicate the extraction procedure, the KK 
relationship between the real and imaginary components of the refractive index must hold at all 
wavelengths. The correct form of Im(N) is easy to obtain from Eq. ([3]i; one simply selects the 
sign such that the imaginary component is positive. We can verify the validity of the similarly 
obtained real components by integrating Eq. (HI over frequencies that effectively span the range 
to °°. In practice, very accurate values for Re(N) may be obtained by numerically integrating 
Eq. (|4| over a very wide, finite range of frequencies. The infinite-frequency limit of the refrac- 
tive index, «oo, is given by the extraction procedure in the limit of very small wavelength. In 
the limit of co ^f°°, all branches in the inverse cosine function converge and all ambiguities are 
removed. 

It should be noted that we have used Kramers-Kronig relations concerning N = n + iK. Some- 
what less convenient for our purposes, but more fundamental, KK relations exist for e and jj.. 
The validity of our KK relations requires that A^ be analytic in the upper half of the complex 
plane, N ^ rioo &s CO —>■ °°, and no conductivity at © = 0. These conditions in relation to metama- 
terials were discussed by Skaar ll30l . In particular, the metamaterial must be a passive medium, 
with no gain. However, it has also been pointed out that A^' has the same analytical properties 
as e and jj. separately, and therefore a KK relation for A^^ would be appropriate in the case of 
gain lISTl . We have verified numerically that the computed n and K from our approach discussed 
below satisfy this other KK relation. In future work, it may be interesting to consider the use of 
this relation directly. 

3. Applications 

We consider a two-dimensional system of gold nanorods separated by Si02 similar to that in- 
vestigated by Shalaev et al. ll20l . The system is illustrated in Fig. Uta), with the parameters 
given in Ref. f20l to obtain negative refractive index; the rods of thickness d=13 nm and width 



w=450 nm are separated by Si02 with a total slab thickness of D=160 nm. The system is peri- 
odic in the x-direction with a periodicity of h=480 nm. We obtain the complex transmission and 
reflection coefficients with finite-difference time-domain (FDTD) simulations as implemented 
in the freely available Meep FDTD package [[27 1. This particular structure, as noted in Ref 1201 
and depicted in Fig [lib), clearly exhibits a negative refractive index in the 850-900 nm region. 
Our present results reproduce well those given in Ref. 1201 . As an aside, we note that the per- 



Fig. 1. The unit cell of a two-dimensional NIM composed of two Au rods separated by 
Si02 (a). The real component of the refractive index (b) is clearly negative in the 850-900 
nm range. 
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meability for this structure is positive; the structure is a NIM but not a double negative material. 
The material does, however, satisfy the necessary conditions given in Ref. ||6l. Optimizing this 
structure for negative real component of the refractive index, we consider first only varying the 
rod thickness, d, while holding all parameters at their values listed above. Figure|2]illustrates the 
real and imaginary components of the refractive index obtained from the extraction in Eq. (O as 
well as transmission and reflection spectra. The real component of the refractive index for the 
d = 10 and d = 14 nm cases (Fig.|2a)), like the d = 13 nm case, is unambiguously determined 
by the extraction procedure. However, a discontinuity develops in the m=Q branch for d = 18 
nm. As stated above, the macroscopic refractive index should be independent of slab thickness, 
and this property should be exploited in assigning branches. In practice, though, a NIM almost 
never exhibits this nice property, and the assignment must be made on the basis of maintaining 
the continuity of the refractive index as a function of frequency. A continuous refractive index 
function can be obtained for d = 18 nm by combining the m=Q and m=l branches, resulting in 
the upper continuous curve. However, this choice gives the refractive index as a discontinuous 
function of the rod thickness, d. One would expect that such an abrupt change in character of 
the refractive index would be associated with an equally abrupt change in either the imaginary 
component of the refractive index or the transmission or reflection spectra, but each of these 
do indeed vary smoothly with the geometric parameter, d. The imaginary component of the 
refractive index does lose a sharp peak when increasing d from 14 to 18 nm, but this change is 
subtle compared to that suggested by Fig.|2la). We propose that the correct form of Re(N) may 
be obtained unambiguously by the integration of Im(N) according to Eq. (|4]i. 

Consider again the well-behaved system (d=13 nm) in which negative refractive index is 
observed in the 850-900 nm range. The values for the index extracted from Eq. ^ and those 
obtained by integration of Eq. (|4]i should be equivalent, provided that we have integrated Im(N) 



Fig. 2. Real (a) and imaginary (b) components of the refractive index extracted according 
to Eq. l3j as well as transmission (c) and reflection (d) spectra for various rod thicknesses, 
d. All other parameters are held constant at the values detailed in the text. 
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over an appropriately large range of frequencies. Table[T]lists the non-parallelity errors (NPE) 
obtained for Re(N) obtained from the KK integration as compared to that obtained from the 
extraction procedure for several different wavelength ranges. NPE are defined here by taking 



Table 1. Non-parallelity error (NPE) in the Kramers-Kronig integration as compared to 
extracted values of Re(N). The NPE is defined as the difference between the most positive 
and most negative errors between the KK and extracted values of Re(N) over the range 
700-1300 nm. 
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the difference between the most positive and most negative errors in the KK integration as 
compared to the extracted values of Re(N) in the range of 700-1300 nm. Integration over only 
the range of interest (700- 1 300 nm) leads to very large errors, and the integrated values of Re(N) 
are virtually meaningless. The KK value for Re(N) become quite accurate when we extend 
the integration range to 100-5000 nm, as is evidenced by the NPE in Table [T] and visually in 
Figs. [3a) and (b). The KK integration reproduces the negative region of refractive index in the 



Fig. 3. Real and imaginary components of the refractive index as obtained from the ex- 
traction in Eq. ([3j and, for the real component, the Kramers-Kronig integration of Eq. l|4j. 
Results are presented for a gold nanorod system with d=13 nm. 
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850-900 nm region as well as all other features from 100-5000 nm. The very large and very 
wide features present in Im(N) in Fig. [2c) illustrate the importance of such a wide range of 
frequencies for obtaining reliable results; these broad features significantly contribute to the 
KK integral. 

Now consider the d=18 nm case in which a discontinuity developed in the m=Q branch, and 
the choice for the correct branch is made subjectively. Forcing Re(N) to be a continuous func- 
tion of wavelength would require one to combine the m=Q and m=l branches to form the upper 
curve in Fig.[2ja), but as shown in Fig.[4l the correct choice as determined by KK integration 
involves only the m=Q branch in the 700-1300 nm range, yielding a negative index with a mini- 
mum around 880 nm. The original branch assignment based solely on the extraction procedure 
and continuity of the refractive index predicts no negative index and is thus qualitatively incor- 
rect. Ambiguous situations such as this where index assignments are based on the continuity of 
Re(N) are surprisingly common, and the KK integration can provide a very easy and necessary 
check as to the validity of branch assignments. The question remains as to why the extraction 
procedure results in a spurious discontinuity in such an important frequency range. In general 



Fig. 4. Real component of the refractive index as obtained from the extraction in Eq. ([3) or 
the Kramers-Rronig integration of Eq. (O for a gold nanorod system with d=18 nm. 
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we observe that, for a gold nanorod NIM of fixed width, D, as the ratio of d to D increases, a 
false discontinuity inevitably develops in the extracted values of Re(N). The KK integration, by 
enforcing causality, can be used to interpolate through these problematic regions and select the 
proper branches on either side. 

The combined extraction/KK procedure for determining refractive indices, in removing am- 
biguities in branch assignments, allows for a more "black-box" approach to the optimization 
of NIM for either very negative refractive index or for a NIM quality factor or figure of merit 
(FOM) related to the ratios of the real and imaginary components if N: 



FOM: 



Re{N) 
' Im{N)' 



(6) 



We have optimized the structure in Fig. [Tta) for the FOM by varying the parameters d, D, and 
w, sampling a total of 768 structures with d, D, and w in the ranges 12-19 nm, 120-230 nm, 
and 400-470 nm, respectively. The structure exhibiting the largest FOM as given by Eq. ^ was 
found to have d=I7 nm, D=I90 nm, and w=430 nm, with a NIM figure of merit, FOM = 1.517. 
This value is nearly double the value of 0.829 that corresponds to the structure suggested in 
Ref. II20I with d=13, D=160, and w=450 nm, and the optimized material is a double negative 
material, where both jx and e are simultaneously negative. Recall that the FOM=0.829 structure 
of Ref. II20I has a positive relative permeability. Fourteen structures were found to have figures 
of merit greater than 1 .4, and, surprisingly, four of these structures might not have been found 
with the simplistic, continuous-function approach to the extraction procedure. Figure |5] depicts 
the m=0 and m=l brances of Re(N) for the optimal FOM=1.517 as well as the best structure 
found that exhibits a false discontinuity in the extracted index, with FOM= 1.442 (d = 17 nm, D 
= 180 nm, w = 440 nm). Clearly, the KK relations allow us to easily and unambiguously select 
the proper branches of Re(N) obtained from the retreival. Amazingly, the high-quality structure 
whose refractive index is represented in Fig.|5jb) would not have been considered a NIM with- 
out the KK analysis. This FOM=1.442 structure, like the optimal FOM=1.517 structure, is also 
a double negative material. 



4. Conclusions 

We have developed a reliable and robust methodology for extracting effective refractive indices 
for metamaterials by combining the usual S-parameter extraction procedures with the concept 



Fig. 5. Real component of the refractive index as obtained from the extraction in Eq. ([3) or 
the Kramers-Kronig integration of Eq. Q for a gold nanorod system with (a) FOM=1.517 
and (b)FOM= 1.442. 
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of causality as expressed in the Kramers-Kronig relationship. Causality absolutely must be 
enforced when extracting effective parameters from the three-layer homogenous metamaterial 
model. Our results suggest that simplistic extraction models in which one copes with branching 
ambiguities by enforcing continuity of the refractive index as a function of frequency may result 
in qualitatively incorrect results. In addition, in regions where Im(N) is very close to zero, the 
sign chosen for Eq. Q may suffer from numerical noise; this sign issue has no implications 
for Im(N), as it is virtually zero, but Re(N) may be given the incorrect sign. Re(N) obtained 
from KK integration will be insensitive to this issue when Im(N) is close to zero and will 
thus yield more reliable results. The Kramers-Kronig relations naturally guide branch selection 
and are simple enough that they may be used in a "black-box" sense for the optimization of 
metamaterial geometries for desired properties. The present methodology was used to optimize 
a gold nanorod system for a NIM quality factor, g, resulting in a structure with a g-f^ctor 
nearly double that of the original starting structure. Several high-quality structures were also 
obtained that would not have been considered NIM with existing extraction procedures. 

The KK relations can also provide a useful gauge as to the validity of the assumption of meta- 
material homogeneity that is central to the S-parameter retrieval procedure. One can view the 
false discontinuities in the refractive index functions as a the first symptoms of the breakdown 
of the validity of the homogeneous three-layer model. In extreme cases, where the amount of 
metal relative to the dielectric present in the sample becomes significant, the KK relation of Eq. 
(|4|i can only qualitatively reproduce the extracted index over the entire frequency spectrum. The 
KK results are not incorrect; the discrepency merely reflects the failure of the assumption of 
the homogeneous nature of the metamaterial. In these rare instances, the extraction procedure 
as a whole may not be an appropriate way to obtain refractive indices. 
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